Anisotropic Neutron Spin Resonance in the Fe-based BaFei 9 Ni .i As 2 Superconductor 
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We use polarized inelastic neutron scattering to show that the neutron spin resonance below T c 
in superconducting BaFe1.9Nio.1As2 (T c = 20 K) is purely magnetic in origin. Our analysis further 
reveals that the resonance peak near 7 meV only occurs for the planar response. This challenges 
the common perception that the spin resonance in the pnictides is an isotropic triplet excited state 
of the singlet Cooper pairs, as our results imply that only the 5*001 = ±1 components of the triplet 
are involved. 

PACS numbers: 74.70.Xa, 78.70.Nx, 75.30.Gw 



I. INTRODUCTION 

The neutron spin resonance is a collective mag- 
netic excitation appearing in the superconducting state 
of high-transition temperature (high-T c ) copper oxide 
superconductors 1 - - —. Since its initial discovery in optimal 
hole-doped YBa 2 Cu306+ a ~ r — , the resonance has been 
found in electron-doped cuprates^, heavy fermion^ 7 -, and 
iron arsenide superconductors^ - — . Below the supercon- 
ducting transition temperature T c , the intensity of the 
resonance increases like the superconducting order pa- 
rameter and its energy scales with Tj±. Although the res- 
onance appears to be a ubiquitous property of unconven- 
tional superconductors 1 - - — , its microscopic origin and re- 
lationship with superconductivity are still debated 1 ^. In 
all these materials, the resonance occurs at the antiferro- 
magnetic (AF) wavevector Q of the parent compound. It 
is thought to be a triplet excitation of the singlet Cooper 
pair o 1 ' 15 , implying a superconducting order parameter 
that satisfies A^+q = — A^. In the iron arsenide super- 
conductors, this condition is satisfied by an order param- 
eter whose sign reverses between the electron and hole 
pockets^ - —. If this picture is correct, one would expect 
that the triplet would be degenerate, and thus direction- 
ally isotropic in space. For the optimal hole-doped high- 
Tc cuprate superconductor YBa 2 Cu306+:r, polarized in- 
elastic neutron scattering experiments suggest that this 
is indeed the case 2 -^. 

We report polarized inelastic neutron scattering results 
for the optimal electron-doped iron arsenide supercon- 
ductor BaFe1.9Nio.1As2 (T c = 20 K; |:: :l . We find that 
the resonance previously observed around 7 meV at the 
AF wavevector Q = (0.5,0.5, 1) (reciprocal lattice units, 
'r.l.u') is entirely magnetic, but displays strong spin- 
space anisotropy, with a peaked response near the reso- 
nance energy occurring only for the planar response. This 



is different from the momentum space anisotropy, where 
the spin correlation length might be different along differ- 
ent crystallographic directions 2 ^ - — . Our results indicate 
a strong spin-orbital/lattice coupling in iron arsenide su- 
perconductors (quite different from the cuprates), and 
are a challenge to the common assumption that the reso- 
nance represents an isotropic singlet-to-triplet excitation. 



II. EXPERIMENTAL DETAILS 
A. Sample 

We chose the iron arsenide superconductor 
BaFe1.9Nio.1As2 because this material has no static 
AF order [Fig. 1(b)], exhibits a well-defined neutron spin 
resonance near 7 meV at Q = (0.5,0.5, 1) above a clear 
spin gap, and is available in large, homogeneous single 
crystal s 10 ' 11 . We co-aligned ~5 grams of single crystals 
(with mosaic of 3° full-width-half-maximum) in the 
(H, H, L) scattering pian o 10 ' 11 , where the wavevector Q 
is indexed Q = Ha* + Kh* + Lc* with a* = a 27r/a, 
etc., a = b = 3.93 A and c = 12.77 A [Fig. 1(a)]. In this 
tetragonal notation, the AF order and resonance occur 
at Q = (0.5, 0.5, L) with L = ±1, ±3, 
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B. Polarized Neutron Analysis 

We carried out polarized inelastic neutron scattering 
experiments using the Cryopad capability of the IN20 
triple-axis spectrometer at the Institut Laue-Langevin, 
Grenoble, France. Neutron polarization analysis is the 
only way to conclusively separate the magnetic sig- 
nal from lattice effects, and to determine the spatial 
anisotropy of the magnetic excitations. 
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FIG. 1: (color online) (a) Crystal structure of BaFe2As2- 
(b) Magnetic and superconducting phase diagram of 
BaFe2-a;Ni :E As2 with the present composition highlighted 
with an arrowi^. (c) Schematic showing a fluctuating atomic 
magnetic moment vector m, and the components m y and m z 
which are probed. Neutron scattering intensity is related to 
the square of the fluctuating moment components in the y and 
z direction, M y = (m y ) and M z = {m 2 z ) respectively (which 
are both defined perpendicular to x, where x is parallel to 
the wavevector Q). (d) Schematic showing crystallographic 
in-plane (1, 1, 0), (f , — f , 0) and out of plane (0, 0, f ) directions 
compared with the x, y, z directions defined above. 



In principle, polarization analysis can be used to com- 
pletely separate magnetic (e.g. spin fluctuation) and nu- 
clear (e.g. phonon) scattering because the spin of the neu- 
tron is always flipped in a magnetic interaction where the 
neutron polarization is parallel to the wavevector transfer 
Q 24 . For convenience, we define the neutron polarization 
directions along Q as x, perpendicular to Q but in the 
scattering plane as y, and perpendicular to Q and the 
scattering plane as z, respectively [Fig. 1(c)]. At a spe- 
cific wavevector and energy, we measured the six cross- 
sections which correspond to the three incoming neutron 
polarization directions x, y, z, with the outgoing neu- 
tron polarization either parallel to the incoming (neutron 
non-spin flip or NSF) or anti-parallel (neutron spin-flip 
or SF) . The measured neutron cross-sections are then ac- 
cordingly written as and tr^ F where a — x,y,z 24: . 
With the Cryopad setup, these cross-sections can be mea- 
sured with the sample in a strictly zero magnetic field 
(< 10 mG), thus avoiding errors due to flux inclusion or 
field expulsion in the superconducting phase of the sam- 
ple. 

We define the magnetic intensity of excitations with 
fluctuating magnetic moments pointing parallel to the 
(1, 1,0) (in-plane) direction as M(110), and the intensity 
of fluctuating moments pointing out of plane as Af(001). 
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FIG. 2: (color online) Energy scans at Q = (0.5, 0.5, 1), show- 
ing crf F (magnetic) and <t£ ,sf (nuclear) scattering for (a) 1.5 K 
and (b) 30 K. (c) (H, H, 1.1) Q-Scan through the resonance 
position, showing <r^ F and er£ ,SF measured at a constant en- 
ergy 7 meV. The narrow peak at (0.625, 0.625) is temperature 
independent spurious scattering, (d) Trajectory in reciprocal 
space of the (H, H, 1.1) scan. Solid lines are guides to the eye 
for all plots except where otherwise stated. 



Our experiment probes M y and M z , the magnetic inten- 
sity of excitations with the moment parallel to y and z 
respectively [see Fig. [UJc)]. Due to tetragonal symmetry 
M (110) = M(110) = M z , and M(001) can be found from 
M y using M y = M(110)sin 2 6> + M(001) cos 2 6>, where 
9 is the angle between wavevectors (1,1,0) and Q [see 
Fig-ETd)]. 

The measured cross-section can be written 
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with a nuclear scattering strength N (containing both 
phonon and inelastic incoherent nuclear scattering), and 
b± and b 2 account for instrumental background (and 
nuclear-spin incoherent scattering). R specifies the qual- 
ity of the neutron beam polarization (so that leakage be- 
tween SF and NSF channels caused by imperfect polar- 
ization are taken into account). In our setup we measured 
R by the leakage of nuclear Bragg peaks into the (mag- 
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FIG. 3: (color online) Energy scans at Q = (0.5,0.5, 1). Raw 
(a) 1.5 K crf F 2 an d (b) 1-5 K 0^«f« cross-section data. Clear 
anisotropy is evident because cr^ F ^ crf F and cr^ SF 7^ cr^ SF . 
(c,d) Raw data taken at 30 K. (e)-(f) In-plane and out of 
plane magnetic response [the extracted M(110) and M(001) 
using raw data in Eq. 2, divided by the Bose factor ('b.f.')] 
at 1.5 K and 30 K respectively. Data in (e) and (f) are also 
corrected for second-order monitor over-counting. 



of neutron polarization direction. To extract M y and 
M(110) from the raw data we can use (from Eq. 1) 
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where c = (R— 1)/(R+ 1). With Eq. 2, we can estimate 
the energy and wavevector dependence of M(110) or M y 
from a weighted average of the pair of values calculated 
using the SF and the NSF data. 



III. RESULTS 

Figs. 2(a)-(b) show crf F (primarily magnetic) and cr^ SF 
(primarily nuclear) energy cuts at (0.5,0.5,1) taken at 
temperatures of 1.5 K (< T c ) and 30 K (> T c ), re- 
spectively. As the temperature decreases, it is clear 
that the nuclear scattering (<x£ ISF ) changes very little 
with temperature, while the magnetic scattering (crf F ) 
around 7 meV is enhanced, and below ^3 meV be- 
comes gappe d 10 ' 11 . These data unambiguously demon- 
strate that the resonance is purely magnetic without any 



lattice contribution. Fig. 2(c) shows a T=1.5 K Q cut 
along the (if, if, 1.1) trajectory at 7 meV [Fig. 2(d)]. 
Consistent with unpolarized measurement s 10 ' 1 1 , the data 
prove that the resonance is magnetic scattering centered 
at (H, K) = (0.5,0.5). 

Having established the magnetic nature of the reso- 
nance, we now probe the anisotropy of the spin fluctu- 
ation spectrum by measuring trf F z and crj! , ?, F and us- 
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ing Eq. 2 to calculate M(110) and M(001). ex- 
clusively probes the in-plane spin fluctuations Af(110) 
and erf F gives the intensity of moments fluctuating along 
My ~ M(001). Finally, of F is the magnetic part of 
the cross-section observed in unpolarized measurements 
and provides the sum of the magnetic scattering, in this 
case My + M(110). For isotropic paramagnetic spin 
fluctuations, one expects M y = Af(110) and this ap- 
pears to be the case for the resonance in optimal doped 
YBa 2 Cu 3 6+ ^. 

Figs. 3(a)-(d) show all six scattering cross-sections 
crf F z and cr^| F raw data taken at Q = (0.5,0.5,1) 
below and above T c . While the resonance at 7 meV 
is clearly seen in the crf F data at 1.5 K [Fig. 3(a)], a 
comparison of <jy and erf F shows that the former has 
a peak while the latter is featureless near the resonance 
energy. Since ~ M(001) and erf oc M(110), these 
data immediately suggest anisotropic spin fluctuations 
near the resonance. By using all six scattering cross- 
sections in Figs. 3(a) and 3(b), we extract the energy 
dependence of Af(110) and Af(001) magnetic scattering, 
and subsequently convert the extracted data to a mag- 
netic response, X110 and X001, [ Fi S- 3(e)] by dividing 
out the Bose population factor [also, we can instead ex- 
tract M(110) and Af(001) from only the three SF cross- 
sections, in which case we get quantitatively very similar 
results]. It is clear that the in-plane response (xi'in) rc ~ 
sembles a peak centered at around 7 meV, whilst the out 
of plane Xqqi nas a much lower energy scale. 

Figs. 3(c) and 3(d) show cr^. F y z and c^y F measured 
at 30 K. Compared with the 1.5 K data, the most obvi- 
ous changes in the erf F z data are the suppression of the 
resonance and the low-energy spin gap. Fig. 3(f) plots 
the energy dependence of the extracted, Bose-factor di- 
vided M(110) and A/(001) at 30 K. In addition to the 
disappearance of the low-temperature spin gap, it can be 
seen that \iio s * m nas a broad peak near E = 7 meV, 
while X001 i s a g am relatively featureless. Comparison of 
the Figs. 3(e) and 3(f) reveals clear evidence for the res- 
onance peak at 7 meV above a spin gap of ~3 meV in 
XiiQj while x'001 ^ s featureless near 7 meV with a spin 
gap of E < 2 meV 10,11 . Previous unpolarized neutron 
scattering measurements found a spin gap value of about 
3 meV at Q = (0.5,0.5, l)^i. Our polarized data are 
consistent with this as well as the unpolarized results^ 
on the same sample if we combine the extracted M (110) 
and M(001) results (See Appendix [B|. 

To further understand the anisotropy of the spin fluc- 
tuations, we carried out constant-energy scans with all 
three crf F , components at E — 2.5, 7, and 11 meV 
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FIG. 4: (color online) (a)-(c) Q-Scans along the (H, H, 1.1) 
direction at 2.5, 7, and 11 meV respectively, with all three 
spin flip cross-sections measured, (d) L-scan at the resonance 
energy. Crosses resemble estimated instrumental background 
points, extracted from the data shown and cr^| F (not shown) 
using Eq. 1 (assuming b\ ~ 62). The solid line shows the 
expected magnetic scattering assuming an Fe 2+ form factor. 



[Figs. 4(a)-(c)]. At 2.5 meV, below the spin gap, there 
is a peak at the in-plane wavevector (0.5,0.5) in crf F ~ 
M(001), whereas for the identical scan ct^f k m(110) 
is featureless. At low Q (H < 0.4) at 2.5 meV, the 
scattering for ct~ and ar F have different backgrounds 
(see Appendix [A}. This problem is not present in the 
energy scans or other Q-scans taken, where the back- 
grounds b\ and &2 must be independent of polarization 
direction. These constant-i? scans are consistent with 
the constant-Q scans in Fig. 3(e), suggesting that the 
spin gap in M(110) is larger than that in M(001). At 
7 meV, there are peaks in both channels at (0.5,0.5), 
but the anisotropy appears to reverse, implying higher 
intensity in the in-plane M(110) direction. Similar data 
are also found for Q-scans at 11 meV [Fig. 4(c)], con- 
sistent with the constant-Q data in Fig. 3. Finally, we 
plot in Fig. 4(d) the L-dependence of the erf F y scattering 
at 7 meV and 30 K. Instead of simply falling off as the 
Fe 2+ magnetic form facto r 25 i 26 , af F peaks near L = 1 
and decreases rapidly with increasing L above the non- 
magnetic background. These results suggest that the res- 
onance in BaFei.gNig !As2 has c-axis modulations similar 
to underdoped BaFe2- a: Ni 2 ;As2^' and is not entirely two- 
dimensional as in BaFei.84Coo.i6As2^. 



IV. CONCLUSION 

We have performed inelastic neutron measurements 
with full neutron polarization analysis to measure the 
magnetic anisotropy of the spin fluctuations in optimally 
doped superconducting BaFe1.9Nio.1As2. We have ob- 
served the magnetic response of the iron spins pointing 
along in-plane [parallel to (1,1,0)] and out of plane [par- 
allel to (0,0,1)] directions to have very different energy 
dependence. For the in-plane response the resonance 
peak was present, whereas the out of plane response was 
reasonably featureless around the resonance energy at 
7 meV. 

Spin-space anisotropy in the zero energy limit has pre- 
viously been reported from NMR data on an underdoped 
hole-doped composition with no magnetic order—, which 
can be explained in terms of the proximity of the compo- 
sition to the ordered parent compound. However, in our 
non-magnetically ordered sample, we see not just low fre- 
quency anisotropy, but a high frequency novel response 
that has different energy dependencies between different 
spin directions. 

The presence of spin-orbital/lattice coupling could ex- 
plain anisotropy in the spin excitations. In pnictides, this 
is reflected in the undoped compound, where the mo- 
ments are locked to the orthorhombic a-axis^r— [along 
(1,1,0) direction in our tetragonal notation]. The ex- 
istence of the resonance solely in the in-plane response 
is a major challenge to the standard theory where the 
resonance is an isotropic triplet excitation of the singlet 
superconducting ground state. To understand the ori- 
gin of our results within the context of this theory^, we 
note that the spin operators Sno and S^q, when acting 
on the spin singlet superconducting ground state, gen- 
erate the £001 = components of the triplet, whereas 
the spin operator £001 generates the £001 = compo- 
nent. Our results therefore imply that the resonance is 
the S001 = ±1 doublet. To understand this microscop- 
ically, we note that in the magnetically ordered phase, 
the Si 10 = component of the triplet would mix with 
the singlet ground state (since the moments are oriented 
along the orthorhombic a axis). In the non-magnetic 
tetragonal state, this would lead to a low energy doublet 
Sno = 0, S110 = 0, which is equivalent to £001 = ±1 (see 
Appendix |C|) . An alternate possibility is that the reso- 
nance is instead a magnon-like excitation that becomes 
undamped because of the opening of the superconduct- 
ing gap 32 ' 33 , though it is not clear to us why this scenario 
would generate a magnetic response that is localized at 
a particular energy. 

The work at UT/ORNL is supported by the U.S. DOE 
BES No. DE-FG02-05ER46202, and by the U.S. DOE, 
Division of Scientific User Facilities. Work at IOP is 
supported by the CAS. Work at ANL is supported by 
the US DOE under Contract No. DE-AC02-06CH11357. 
OJL and TE were supported by the DOE BES EPSCoR 
Grant DE-FG02-08ER46528. 
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Appendix A: Neutron Polarization Independent 
Backgrounds and the 2.5 meV Q-cut 

As implied by Eq. 1 of the paper, in principle the back- 
ground scattering into the detector should be the same 
with neutron polarization in x, y, z for any given SF (or 
NSF) measurements since the axes of the instrument do 
not move. However, there is a moving part that does 
change with neutron polarization direction, and that is 
the 'dipole magnet' in the outgoing beam, which rotates 
around the scattered beam axis (with a position depend- 
ing on polarization direction, as well as Q and E) and 
creates the neutron guide field that defines the neutron 
polarization direction. The problem occurs when a choice 
of Q and E conspires to cause both a scattering angle 
that is small, and a dipole magnet position close to the 
horizontal for a certain neutron polarization direction. 
Neutrons can then scatter in grazing incidence from the 
dipole magnet shielding, which can increase the back- 
ground in the detector for that configuration over other 
neutron polarization directions. 

At low Q (H < 0.4) at E =2.5 meV, these problematic 
conditions are satisfied creating an extra background in 
the cr£ ISF,SF configurations. However, the dipole magnet 
is away from horizontal at the same energy and wavevec- 
tor for the cr^ SF:SF and cr^ SF,SF configurations, which 
therefore have lower backgrounds. At higher energies 
(and other wavevectors), the dipole magnet is never close 
to horizontal when the scattering angle is small enough 
for a grazing incidence to reach the detector, and so there 
is no difference between backgrounds of different neutron 
polarization configurations. We have confirmed that this 
is indeed the case, by comparing backgrounds extracted 
for all the data collected, and found an anomalous effect 
only for the low Q region at 2.5 meV. 

In conclusion, at H < 0.4 in the 2.5 meV Q-scan there 
may be a difference between backgrounds in configura- 
tions with different neutron polarizations (and thus, in 
this case the assumption in Eq. 1 may not be valid). 
However, this is not a problem in any other scans, and 
most importantly does not affect the energy scans at any 
point. Therefore, as required to correctly extract M(110) 
and M(001), the assumption that the background is neu- 
tron polarization direction independent is a good one for 
the energy scans. 



Appendix B: Comparison of Extracted Data and 
Previous Unpolarized Results 

From the present study, from observing the two differ- 
ent spin gaps at 3 meV and < 2 meV, and different max- 
ima at approximately 7 and 3 meV in the M(110) and 
M(001) channels, one might expect to see these feautures 
in unpolarized data. The same compound has been previ- 
ously studied^! by unpolarized neutrons in the (H, K, 0) 
scattering plane (different from the scattering plane in 
the present study). Although there is a resonance at 
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FIG. 5: (color online) Intensity expected for energy scans at 
(0.5, 0.5, 0) for an unpolarized experiment, calculated using 
M(110) and M(001) from the present polarized data. 



7 meV and a spin gap around 3 meV, the dynamic sus- 
ceptibility does not have a peak near 3 meV. Here we 
show that these results are entirely consistent with the 
present polarized neutron scattering results. 

In the unpolarized experiment, the magnetic scatter- 
ing measured at (0.5, 0.5, 0) is proportional to M(110) + 
Af(001) for the crystal alignment used. If we assume 
minimal L dispersion, then we can take the M(110) and 
M (001) values from our present study (where L = 1) and 
simulate the (L = 0) unpolarized data with no unknown 
parameters. We can then compare our simulation with 
the experimental data from the unpolarized experiment. 
As can be seen from Fig. [51 the low energy features in 
Af(001) near 3 meV do not cause low energy features in 
the total unpolarized intensity M(110) + Af(001). The 
resulting form of Fig. [5] is consistent with the data in 
unpolarized measurements (in Ref. H3), though the res- 
olution in the unpolarized experiment was much better, 
leading to a much sharper resonance in that study. 



Appendix C: Origin of the doublet resonance 

The spin singlet Cooper pair wavefunction is a product 
of states of the form 



I**, >= \k f, -k i> -\k |, -k t> 



(CI) 



We operate on this state with the spin operator, S, which 
is the sum of S\ and S2 where 1 and 2 denote the two 
electrons of the pair. For the spin raising operator, we 
find 

S+(q)\9 k >= |fc t, -£ + <7 t> ~\k + q t, ~k f> (C2) 

This is the S z =l component of a triplet pair with center 
of mass momentum q (the minus sign being a reflection 
of fermion antisymmetry). Similarly, SL generates the 
S z =-1 component. Had we operated with S z instead, we 
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would have obtained the 5*2=0 component of the triplet. 
Therefore, for a quantization axis along c, Xaa an d Xbb 
generate the 5 C =±1 doublet, whereas Xcc generates the 
S c =0 state. Since we find no resonance response for Xcc, 
the resonance is the S c =±l doublet. To better appreciate 
this result, assume that superconductivity and antifcrro- 
magnetism coexist, corresponding to the spin resonance 
being at zero energy. If one pairs electrons using anti- 
ferromagnetic eigenstates, and then rewrites these pairs 
in terms of paramagnetic eigenstates, the resulting pair 
state is well known to be a mixture of a singlet and the 
S z =0 component of a triplet^ with z parallel to the di- 



rection of the Neel vector. In the isotropic case, the Neel 
vector can point in any direction, which is why the reso- 
nance is a triplet. But for the antiferromagnetic ground 
state of the pnictides, the spins are locked to the or- 
thorhombic a axis. Therefore, the mixed triplet compo- 
nent of the pairs for a coexisting state would be S a =0. If 
we then average in the plane so as to restore tetragonal 
symmetry, then the Sb=0 component would be involved 
as well. Thus we obtain a doublet. If we now rotate the 
quantization axis to be along the c direction, this doublet 
corresponds to 5 C =±1. 
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